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We proof the existence of two different topological classes of low-energy k ■ p Hamiltonians for 
Chern insulators. Using the paradigmatic example of single-valley two-band models, we show that 
k ■ p Hamiltonians that we dub local have a topological invariant corresponding precisely to the Hall 
conductivity and linearly dispersing chiral midgap edge states at the expansion point. Non-local\i ■ p 
Hamiltonians have a topological invariant that is twice the Hall conductivity of the system. This 
class is characterized by a non-local bulk-edge correspondence with midgap edge states appearing 
away from the high-symmetry k ■ p expansion point. 
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Introduction - The discovery of topological insula¬ 
tors (TIs) [Tl[2] has brought to light a new state of quan¬ 
tum matter that has had a tremendous impact in the field 
of fundamental condensed matter physics as well as for 
potential application in spintronics and quantum com¬ 
putation. The TIs are insulating in the bulk, but they 
do possess metallic edge states w, which are topolog¬ 
ically protected against generic perturbations preserving 
the symmetries of the underlying topological class [6] and 
the intrinsic insulating behavior. Most importantly, the 
robustness of these edge states is encoded in a topological 
invariant classifying the ground state of the system: this 
is the celebrated bulk-edge correspondence [3 [S] . With 
non-interacting lattice Hamiltonians at work, the band 
structure of an insulator can be viewed as a mapping 
from the periodic Brillouin zone, which has the topol¬ 
ogy of a torus T‘^ in d dimensions, to the space of Bloch 
Hamiltonians with an energy gap. In two-dimensional 
(2D) systems, this mapping, in turn, allows for the defini¬ 
tion of a topological invariant, which is referred to as the 
Chern number C, without invoking additional symme¬ 
tries such as time-reversal [3] or point-group symmetries 
of the lattice [3]. 

This becomes manifest in a simple fashion by partic¬ 
ularizing to the situation of a two-band Hamiltonian, 
which can be represented in the space of Pauli matri¬ 
ces T as PL (k) = e(k)l 2 — d(k) • r, where, since e(k) 
multiplies the identity matrix I 2 , the eigenstate struc¬ 
ture, and consequently the topological character of the 
electronic ground state of the system, depends only on 
d = {dx, dy, dz). The Chern number of the Hamiltonian 
above is given by 

C = ^ J dkxdky d ■ {dk,,d X dkyd) (1) 

where d is the normalized d vector of unit length, which 
lives on a sphere and thus defines a mapping from the 2- 
torus of the Brillouin zone to the 2-sphere S^. With 
this, C counts the number of times the image of the map¬ 
ping wraps around the sphere, which is obviously an in¬ 


teger, cannot change under smooth deformations of the 
mapping, and thus defines a topological invariant of the 
system. 

Instead of the full electronic band structure over the 
entire BZ one often uses the effective, long-wavelength 
description of a system: its k • p Hamiltonian [10] that 
corresponds to the continuum description of its low- 
energy bands. The obvious question that arises is 
whether and how the topological invariant classifying a 
Chern insulator can be defined in such an effective con¬ 
tinuum description of the system. This is indeed partic¬ 
ularly timely in view of the intense research effort on 
closely-related topological states of matter in inverted 
semiconductors, such as HgTe quantum wells El , or 
bismuth-based materials m- On top of this, effective 
k • p theories have been shown to be essential in order 
to gain insights into the fundamental electronic proper¬ 
ties of the topologically protected surface or edge states 
[IMS]. In this Letter, we address this question and 
provide a classification of low-energy, long-wavelength 
Hamiltonian for 2D Chern insulators, thereby identify¬ 
ing two different topological classes defined in terms of 
topological invariants different in nature. The local k • p 
topological class corresponds to the known case of low- 
energy continuum theories with the momentum space 
one-point compactified to a 2-sphere S^. The topolog¬ 
ical information is then contained in a 5^ —^ mapping 
- with the topological invariant equal to the Chern num¬ 
ber - and manifested in the ubiquitous appearance of 
midgap linearly dispersing chiral edge states with left- 
movers and right-movers crossing each other precisely at 
the k • p expansion point. In the second non-local k • p 
class, the momentum manifold is instead holomorphic to 
a non-compact hemisphere. While this feature has previ¬ 
ously been associated with marginal topological behavior 
m and a consequent absence of topological well-defined 
quantities [muHi, we show that a one-dimensional topo¬ 
logically non-trivial mapping from a circle to a circle is 
still allowed. The ensuing topological invariant then cor¬ 
responds to the topological charge of the equator vortex 
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FIG. 1: (Color online) Top view patterns of the d vector con¬ 
figurations on a Bloch sphere( local k ■ p topological class (a)) 
and a Bloch hemisphere {non-local'k ■ p topological class (b)). 
The two topological classes are dehned in terms of different 
topological invariants. 
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FIG. 2: (Color online) (a) Energy bands for the BHZ tight- 
binding model of Ref. |Tl] in a ribbon of width W = 200a 
with a the square lattice constant. The thick lines represent 
the analytical results for the dispersion of the topologically 
protected edge states as obtained from the low-energy k • p 
expansion at the centre of the BZ. (b) Enlargement in the 
bulk gap region showing a perfect agreement between the an¬ 
alytical and the numerical results. 


and amounts to twice the Hall conductivity of the sys¬ 
tem. An explicit solution of the low-energy theory in 
semi-infinite systems shows, independent of the bound¬ 
ary conditions, the appearance of edge modes at the top 
(bottom) of the valence (conduction) band and vanishing 
Fermi velocity at the k • p expansion point. The non¬ 
trivial topology then yields a non-local bulk-edge corre¬ 
spondence: although the bulk topological invariant can 
be read off from the properties of low-energy electrons, 
chiral edge states of finite Fermi velocity connecting the 
valence to the conduction band are encountered at BZ 
points away from the k • p expansion point. 

Local k • p topological class - When dealing with a 
2D long-wavelength Hamiltonian, the momentum man¬ 
ifold defining the band structure of an insulator corre¬ 
sponds to the real plane Ti?, which is known to be an 
element of the trivial class [ 19 ] - mappings between the 
non-compact and the compact 2-sphere are gener¬ 
ally topologically trivial m . A loophole for this state¬ 
ment appears whenever the infinite k plane can be com- 
pactified as a Riemann sphere m- Put in simple terms, 
this occurs if the normalized d vector converges in the 
k —>■ oo limit to a single point of S^, say the south 
pole, regardless of the direction along which infinity is 
approached, i.e. lim|i(.|_>oo d = const. Mappings from 
this one-point compactified k plane to the 2-sphere sub¬ 
tended by d are topologically well-defined, thereby allow¬ 
ing for a safe definition of the Chern number. To show 
this, we consider the example of the Bernevig-Hughes- 
Zhang (BHZ) model [TT] originally introduced to discuss 
the time-reversal Quantum Spin Hall insulating phase in 
HgTe/CdTe quantum wells. Assuming a 17 ( 1 ) spin sym¬ 
metry, each of the two spin sectors of the BHZ Hamilto¬ 
nian corresponds to a long-wavelength k • p theory of a 
Chern insulator with the d vector reading 

( 2 ) 


where Sz = ±1 distinguishes the two spin sectors. For 
m > 0 the normalized d vector visits both the south and 
the north pole, wraps the unit sphere once, thereby im¬ 
plying an integer Chern number C = Sz- In the opposite 
m < 0 regime instead, the image of the mapping does not 
cover the full sphere yielding a trivial (7 = 0 Chern num¬ 
ber. The well-defined topological properties of this long- 
wavelength Hamiltonian are reflected in the ubiquitous 
emergence of topologically protected edge states within 
the k • p framework. To show this, we consider the long- 
wavelength model Hamiltonian defined on the half-plane 
given by a: > 0 and consider boundary conditions that 
ensure a self-adjoint extension of the half-plane Hamilto¬ 
nian, i.e. = 0 with the domains of LL 

and T-0 equal. Those correspond to the usual Dirichlet 
fixed boundary conditions (FBC) [20] i which clamp the 
wave function to zero at the boundaries, and to ” natural” 
boundary conditions (NBC), that can be also derived by 
minimizing a variational energy functional [21] , explicitly 
reading dH/dkz^iK -t -idx)'ip\x=o = 0- 

We then look for topologically protected edge states 
by assuming the ansatz for the localized wavefunction 
^ 9I(Ai) < 0. By solving the 

corresponding Schrodinger equation, we obtain midgap 
chiral edge states of opposite chirality in the two spin sec¬ 
tors with energy E = —Szky. The corresponding wave- 
functions read 

^oc[i,-Szf ( 3 ) 

where the ± are for the NBC and the FBC boundary 
conditions respectively and the Ai _2 satisfy 

Ai,2 = -^±y^+fcy-m, (4) 

which imply a renormalizability of the edge states wave- 
function in the bulk gap region m > ky. The dispersion 
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of the edge states is shown with the thick lines in Fig. 
We find an excellent agreement with the numerical re¬ 
sults obtained by solving the square lattice tight-binding 
Hamiltonian of Ref. HH in a ribbon geometry, explicitly 
proving that the absence of a short-distance cutoff in the 
long-wavelength Hamiltonian does not change neither the 
topological properties nor the electronic characteristics of 
the topologically protected edge states. 

Non-local k • p topological class - Asa matter of fact, 
for certain Chern insulating lattice models a one-point 
compactification of the effective low-energy k • p theory 
can strongly modify the intrinsic topological properties. 
This can be explicitly shown by considering a lattice 
model for spinless fermions on the checkerboard lattice 
[22] in the quantum anomalous Hall (QAH) state |23|. 
The corresponding d vector of the lattice Hamiltonian 
can be written as [21| 

( 4t cos (j) cos ^ cos ^ \ 

At sin (j) sin ^ sin ^ 1 (5) 

2 t 2 (cos kx — cos ky) J 

where t and t 2 are the nearest-neighbor and next-nearest- 
neighbor hopping amplitudes, while ±(/> are local mag¬ 
netic fluxes picked up by the electrons while hopping 
around each plaquette clock and aniclockwise respec¬ 
tively. Without loss of generality, we will assume in the 
following t,t 2 >0 and 0 < (j) < tt/2. The effective low- 
energy, long-wavelength Hamiltonian can be derived ex¬ 
panding Eq. close to the corner of the Brillouin zone 
M = {tt, tt} as 

d ~ {tcos(j)kxky,N{'k),t2{kl - k^)), (6) 

where the time-reversal symmetry breaking mass A(k) = 
Aq — Aik^. Due to the absence of a compact k-space, 
the mapping d is topologically trivial, reflected in the 
fact that the ensuing Chern number takes a meaningless 
non-integer value as long as Ai ^ 0 . 

We therefore first seek for a natural one-point com¬ 
pactification by assuming a modified functional form of 
the time-reversal symmetry breaking mass A(k) = Aq — 
Aik^ — A 2 k'^. Since in this case d defines a topological 
non-trivial mapping, integers Chern numbers have to be 
restored. And indeed we find C = Ao/|Ao| -I- A 2 /IA 2 I G 
{± 2 , 0 } as easily follows from the fact that the d-wave 
symmetry of the effective low-energy Hamiltonian implies 
that the unit sphere subtended by d can only be wrapped 
twice. This result is in striking contrast with the analysis 
of the checkerboard lattice model Eq. , which predicts 
C = ±1, even though the integrand in Eq. j^is strongly 
localized close to the M point of the Brillouin zone [c.f. 

Figi- 

Having established that the topological characteristics 
of the checkerboard lattice model cannot be correctly de¬ 
scribed within a one-point compactified k • p Hamilto¬ 
nian, we now proceed to introduce the second topolog¬ 
ical class of effective long-wavelength Hamiltonians for 



FIG. 3: (Color online) Plot of the integrand of Eq.[^as a func¬ 
tion of momentum k over the entire Brillouin zone, obtained 
from the checkerboard lattice model in the quantum anoma¬ 
lous Hall state with the parameter t 2 /t =1/2 and </ = 1/10. 
The integrand peaks close to the M points. 


Chern insulators and the topological invariant that can 
be defined thereof. We start out by requiring the momen¬ 
tum manifold G to be holomorphic to a non-compact 
hemisphere. This can be accomplished if in the k —^ 00 
limit the normalized d vector is planar and thus con¬ 
fined to the equator of a 2-sphere . It then follows 
that for a closed loop in momentum space in the large 
|k| limit, the planar d vector defines a topologically non¬ 
trivial mapping from a circle 5^ to a circle 5^, whose 
ensuing topological invariant counts the number of times 
d encircles the origin. 

In order to show how this winding number is related to 
a well defined physical response, we recall that the Hall 
conductivity an in units of e^//i can be written as a line 
integral over a closed loop C in momentum space 

^ £dk- A{k), ( 7 ) 

where the Berry connection of the occupied lower band 
A(k) = —i (M(k)l V |u(k)) in terms of the Bloch wave- 
functions |u(k)) |25|. Since the Berry phase in Eq. 
is gauge dependent, however, the gauge needs to be 
fixed. We accomplish this by requiring that the wave- 
functions at the dz = 0 equator are able to contract into 
the north or south pole. With this, it follows that the 
lower band eigenstate at the dz = 0 equator has to be 
written as ■u(k) = {l,dx + idy)'^ 1^/2 for ^3 > 0 and 
M(k) = {doc — idyA )'^for ^3 < 0. This gauge fix¬ 
ing procedure, in turns, yields an = Wdz/{2\dz\), where 
W = l/(27r) dk ■ [do;Ndy — is a topological in¬ 

variant corresponding to the topological charge of the d 
vortex [ 2 H] . 

For the effective low-energy theory of the checkerboard 
lattice model Eq. the momentum manifold can be ren¬ 
dered to be holomorphic to a hemisphere by neglecting 
the momentum dependence of the time-reversal symme¬ 
try breaking mass A(k) = Aq. With this, it follows that 
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FIG. 4: (Color online) (a) A checkerboard lattice (b) Energy 
bands for ribbons in the QAH state with = 1.2, cj) = 
1/10. The edge state crossing always occurs away from the 
low-energy k ■ p expansion point M=(7r, tt) where the bulk 
band gap is smallest, (c) Enlargement close to the M point 
projection. The dashed line are the edge states as obtained 
from the low-energy k • p theory, (d) Behaviour of the chiral 
edge states Fermi velocity vanishing at the k ■ p expansion 
point. 


the Hall conductivity an = Ao/|Ao|, which is in per¬ 
fect agreement with analysis of the lattice Hamiltonian, 
and reconciles the topological properties of the long- 
wavelength theory with the short-distance cutoff one. 

Next, we prove that the two topological classes of con¬ 
tinuum k • p Hamiltonians can be discriminated not only 
from their different topological invariants but also from 
the ensuing properties of the topologically protected edge 
states. This can be immediately seen by analyzing as 
before the long-wavelength Hamiltonian of Eq. on a 
half-plane with FBC and NBC boundary conditions and 
scanning for localized wavefunctions exponentially decay¬ 
ing into the bulk. It can be then shown that such edge 
states are non-dispersive and located at the top (bottom) 
of the valence (conduction) band E = Aosign(fcj,) with 
the wave functions reading 

(x[i,sign{ky)f (8) 

where, similarly to the BHZ model, the ± are for the 
NBC and the FBC boundary conditions respectively and 
the Ai ^2 satisfy 


Ai,2 — — 


tcoscj) 


\ky \ i ky\ 


1 COS^ (j) 


- 1 . 


(9) 


Fig. I^b) shows the energy bands in the QAH state for 
a checkerboard lattice ribbon. At the the low-energy M 
point projection, we find edge states precisely at the top 
(bottom) of the conduction band [c.f. Fig. Qc)] of van¬ 
ishing Fermi velocity [c.f. Fig. |^d)], in perfect agree¬ 
ment with the foregoing k • p analysis . Most impor¬ 
tantly, these edge states acquire a linear dispersion mov¬ 
ing away from the low-energy k • p expansion point and 
eventually connect the valence to the conduction band, as 
guaranteed by the bulk-edge correspondence. This fea¬ 
ture shows the most peculiar characteristic of the non¬ 
local k • p topological class: even though the topological 
invariant can be inferred from the low-energy bulk band- 
structure, it exists a non-local bulk-edge correspondence 
predicting midgap edge states at different BZ points. 
This is highlighted by the fact that right- and left-moving 
modes cross each other at an high-symmetry point of 
the one-dimensional BZ different from the low-energy BZ 
point. 

Conclusions - To sum up, we have identified two dif¬ 
ferent topological classes for single-valley k • p Hamilto¬ 
nians of Chern insulators. In the first local k • p topologi¬ 
cal class, the momentum manifold can be one-point com- 
pactified to a Riemannian sphere, and allows for the usual 
definition of the Chern invariant. In the second non-local 
k • p topological class, instead, the momentum manifold 
is holomorphic to a non-compact hemisphere but still 
allows for a one-dimensional topological invariant con¬ 
tained in a mapping from a circle to a circle. These two 
topological classes are characterized by a bulk-edge cor¬ 
respondence different in nature. While in the local k • p 
class the topologically protected edge states appears in 
close proximity to the low-energy points of the BZ, in 
the non-local k • p class the bulk-edge correspondence is 
intrinsically non-local: midgap edge states appear away 
from the k • p expansion point. 
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